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Abstract: Contrary to the customary thought prevailing for long, 
the time reversibility associated with beam-to-beam collisions does 
not really exist. Related facts and consequences are presented. The 
discussion, though involving simple mathematics and physics only, 
is well-related to the foundation of statistical theory. 

1 Introduction 

More than one hundred years ago, the debate concerning time reversibil- 
ity arose in a confusing way: Boltzmann derived his kinetic equation from 
the time reversibility of mechanics while the equation itself was of time 
irreversibility [1]. Even today, though a long time has passed and count- 
less papers in the literature have revealed a variety of aspects of the issue, 
paradoxical things still bother some of us [2]. 

Here, it will be shown that the real problem of Boltzmann's theory is re- 
lated not to the time irreversibility assumed by it, but to the time reversibility 
assumed by it. To make the topic more intriguing and more profound, the in- 
vestigation will manifest that any attempt to formulate the 'true distribution 
function' will fail in the ultimate sense. 

2 Particle-to-particle and beam-to-beam col- 
lisions 

Before entering the detailed discussion, it is essential to establish distinction 
between particle-to-particle collisions and beam-to-beam collisions. 

A particle-to-particle collision involves two individual particles. The time 
reversal symmetry of it has been fully elucidated in classical mechanics and 
it says that if the collision expressed by (vi,V2) — * (v^, v^) is physically 
possible, where v 1; v 2 are respectively the velocities of the two particles 
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before the collision and v^, \' 2 after the collision, then the inverse collision 
expressed by (— v^, — V2) — > (— Vi, — v 2 ) is also physically possible. This kind 
of time reversibility is not truly relevant to the subject herein and we shall 
not discuss it too much in this paper. 

Beam-to-beam collisions, involving particle densities or distribution func- 
tions, are of great significance to statistical mechanics. For instance, Boltz- 
mann's theory treats f(\)dv as a beam and derives its formulation on the 
premise that certain types of time-reversibility are there. 

However, we happen to realize that no time-reversibility of any form can 
be defined in the context of Boltzmann's theory. This conclusion is surprising, 
seems very imprudent and directly contradicts what has been embedded in 
our mind. In view of such strong resistance, it is felt that a very clear and 
very detailed discussion should be given. In this section the subject will be 
studied intuitively and in the next two sections mathematical investigations 
will be presented. 



Figure 1: A candidate for time reversibility of beam-to-beam colli- 
sion: (a) the original collisions; and (b) the inverse collisions. 

Take a look at Fig. 1. Fig. la shows a process in which two beams with 
two definite velocities collide and the produced particles diverge in space. 
Fig. lb illustrates the inverse process, in which converging beams collide 
and the produced particles form two definite beams. In no need of detailed 
discussion, we surely know that the first process makes sense in statistical 
mechanics while the second one does not. 




(a) 



(b) 
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3 No time-reversibility in terms of cross sec- 
tions 



In the textbook treatment, the time reversibility concerning beam-to-beam 
collisions is expressed in terms of cross sections [1]: 

cr(vi, v 2 -> vi, v 2 ) = <r(vi, v 2 -> V!, v 2 ), (1) 

where <t(vi, v 2 — > v^, v 2 ) is defined in such a way that 

TV = fj(vi, v 2 ->■ v' 1; v 2 )rfvirfv 2 (2) 

represents the number of type-1 particles emerging, after collisions, between 
v^ and v[ + dv[ per unit incident flux and unit time, while the type-2 particle 
emerges between v 2 and v 2 + dv' 2 . 




(a) (b) 



Figure 2: Constraints imposed by the energy and momentum con- 
servation laws. 

Unfortunately, the cross section in (1) and (2) is ill-defined. Notice that 
the energy and momentum conservation laws state that (assuming every 
particle to have the same mass for simplicity) 

c = c' and |u| = |u'| = u, (3) 

where 2c = vi + v 2 , 2c' = v' x + v 2 , 2u = v 2 — vi and 2u' = v 2 — v^. 
Fig. 2a shows how vi and v 2 determine c and u, while Fig. 2b shows how 
c and u = |u| form four constraint constants on v[ and v 2 . Notably, v^, 
as well as v 2 , falls on the spherical shell S of radius u in the velocity space, 
which will be called the accessible shell. With these constraints in mind, two 
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problems associated with the definition (2) will surface by themselves. The 
first is that, after dv[ is specified, specifying dv' 2 in the definition is a work 
overdone (since v[ and v' 2 are not independent of each other). The second 
is that the cross section should be defined in reference to surface elements 
rather than to volume elements. 

To see the second problem aforementioned more vividly, let's consider 
dv[ shown in Fig. 3a, which is cube-shaped with equal sides I. If we let p 
denote the area density of particles on the accessible shell caused by unit 
flux of type-1 particles, the number of type-1 particles found in dv[ can be 
expressed as N « pi 2 . Then, the cross section defined by (2) is equal to, with 
dV 2 omitted, 

N pP_ p 
~ " 7' 



a = 



(4) 



dv[ P 

which depends on / and tends to infinity as the cube becomes smaller and 
smaller. Nevertheless, if d\[ is chosen to be a slim box in Fig. 3b and the 
box becomes slimmer and slimmer, then a tends to zero; if dw\ is chosen to 
be a short box in Fig. 3c and the box becomes shorter and shorter, a tends 
to infinity again. These representative examples inform us that the cross 
section defined by (2), and thus the time reversibility expressed by (1), does 
not mean anything. 



dw[ 



(a) 



(b) 



(c) 



Figure 3: In defining the cross section, the velocity volume elements 
dv[ may take on different shapes. 

As a matter of fact, the above problem can be examined more briefly. For 
any definite velocities vi and v 2 , the six components of v' x and \' 2 are under 
four constraint equations imposed by the energy-momentum conservation 
laws. This literally means that the free space of v' x and is just two- 
dimensional and there are no enough degrees of freedom allowing us to define 
a cross section in reference to six-dimensional volume elements. 
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4 No time-reversibility in terms of velocity 
volumes 



Another form of time reversibility is simultaneously employed in such textbooks[l]: 

dvxdv-2 = dv'-ydv^. (5) 

It should be noted that there is a conceptual conflict between (1) and (5). 
In connection with (1), when incident particles have two definite velocities 
the velocities of scattered particles are allowed to distribute over almost the 
entire velocity space; whereas, in connection with (5), an infinitesimal ve- 
locity range of incident particles strictly corresponds to another infinitesimal 
velocity range of scattered particles. Nevertheless, for purposes of this paper, 
we shall leave this conflict alone. 

The mathematical proof of (5) goes as follows. First, 

dv 1 dv 2 = || J\\dcdu with ||J|| = 

Then, 

dv^dv^ = || J'\\dc'du with \\J'\\ = 
In view of that 

||J|| = ||J'||, c^c' du = du, (8) 

we obtain (5). 

Unfortunately again, the formulation given above also involves errors. All 
equations in (6), (7) and (8) hold except du = du'. Referring to Fig. 2, we 
find that when u is a definite vector, u' distributes over a spherical shell, 
pointing in any direction. This means that if du = u 2 dudfl u is an infinitely 
thin line-shaped volume element (say, dfl u is infinitesimal while du finite), 
the corresponding volume element du' will be a spherical shell with finite 
thickness. It is then obvious that the two elements are not equal in volume. 

The issue can be analyzed more economically in terms of variable trans- 
formation. If we identify vi x ,vi y ,vi z and V2 X , V2 y , viz as six variables and 
identify v' lx ,v' ly ,v' lz and v 2x ,V2y,v 2z as six new variables, then we have 

dvidv 2 = ||J||dvidV 2 . (9) 



<9(vi, v 2 ) 



<9(c,u) 



(6) 



9(vi,v 2 ) 



d(c', u') 



(7) 
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What has been proven by equations (6), (7) and (8) is nothing but || J\\ = 1. 
However, in order for (9) to make sense, there must exist six independent 
equations connecting those variables. In our case, we have four equations 
only. That is to say, the 'variable transformation' is incomplete and expres- 
sion (9) is not truly legitimate. 

It is instructive to look at the issues discussed in the last and this sec- 
tions in a unifying way. If there were no single constraint equation, defining 
the cross section could make sense; if there were six independent constraint 
equations, defining the Jacobian would be meaningful. Since there are four 
and only four constraint equations, neither the cross section nor the Jacobian 
can be defined. 

5 Formulation of beam- to- beam collisions 

It is now rather clear that the concepts and methodologies of Boltzmann's 
theory are in need of reconsideration. 

For instance, one of the major steps in deriving Boltzmann's equation is to 
identify f(v' 1 )dv[ and /(v^rfv^ as two definite beams and then to determine 
how many beam-1 particles will emerge between Vi and vi + dv\ due to 
collisions of the two beams. As has been shown, this context leads to nothing 
but an absurd result: the number of such emerging particles actually depends 
on the size and shape of dv±, varying drastically from zero to infinity. 

In what follows, we shall propose a new context to do the job. Surpris- 
ingly, the formulation will reveal some of deep-rooted properties of statistical 
mechanics. 

To involve less details, we adopt the following assumptions: (i) The 
zeroth-order, collisionless, distribution function of the gas is completely known 
(ii) Each particles, though belonging to the same species, is still distinguish- 
able (which is possible in terms of classical mechanics), (iii) No particle 
collides twice or more. (General treatments can be accomplished along this 
line [3].) 

Referring to Fig. 4, we consider that two typical beams, denoted by 
/i°^( v i)^ v i an d fi I ( v 2)^ v 2) collide with each other, and suppose that a 
particle detector has been placed somewhere in the region. Let AS be the 
entry area of the detector and ANi be the number of the beam-1 particles 
entering the detector within the velocity range Av \v\ Afii during dt. Since 
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any beam of the system can be regarded as the first beam, or the second 
beam, aforementioned, the total distribution function due to collisions is, at 
the detector entry, 

/«( f ,r,v„AS,A„„ A*) „ m 

in which r is the representative position of AS* and vi is the representative 
velocity of Ai^i^Af^. According to the customary thought, when AS, Avi 
and Af^ shrink to zero simultaneously, expression (10) stands for the 'true 
distribution function' there; for reasons to be clear a bit later, we shall, in 
the following formulation, assume that AS and Avi are infinitely small while 
AQi is kept fixed and finite (though rather small). The spatial region — Af^, 
that has been shaded in Fig. 4 and 'opposite' to the velocity solid-angle range 
Afii in (10), will be called the effective cone. It is intuitively obvious that 
the particles that collide somewhere in the effective cone and move, after the 
collision, toward the detector entry along their free trajectories will contribute 
to ANi. (Even if such particles are allowed to collide again, some of them 
will still arrive at the detector freely, which means the concept of effective 
cone holds its significance rather generally.) 



-AQ 1 




I | Detector 



Figure 4: Two beams collide and a particle detector is placed in 
the region. 

Observing the colliding beams in the center-of-mass reference frame, we 
find that the number of collisions in a volume element dr', which is located 
inside the effective cone — Afl, can be represented by, as in Boltzmann's 
theory, 

[rfr'^KJrfv'Jt/fK)^]^^^,^^^], (11) 

where Q c is the solid angle between u' and u, and <7 c (u', u) is the cross section 
in the center-of-mass frame. By integrating (27) over the effective cone and 
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taking account of all the particles that are registered by the detector, the 
right side of (10) becomes 



where Af2 is the solid-angle range formed by a representative point in dr' 
(as the apex) and the detector entry area AS (as the base). In view of 
that AS is truly small and AQi is fixed and finite by our assumption, we 
know that AQq <C AQi and every particle starting its free journey from 
the effective cone and entering the detector can be treated as one emerging 
within Af^. Then, we can rewrite expression (12) as, with help of the 
variable transformation from (v'^vQ to (c',u'), 

,y\ «./ rf cW rf ui l jr» u)/!0, ( c '- u '^ , ( c ' + u ') 



-AQi J AviAQo J J 1 (| r - T , \ 2 ASl v 1 )(vfAv 1 ASl 1 ) 

(13) 

in which ||J|| = d(v[, v' 2 )/d(c', u') and the subindex Av 1 Afi 1 has been re- 
placed by AviAQ . In view of the energy-momentum conservation laws, we 
arrive at 

/ dr' f dc' f dVl u , f dVt c f du 

~ , 2ua c (u>, u)A (0) (c' - u')/ 2 (0) (c' + u') (14) 



■J 



(| r _ r '\ 2 An v 1 )(vlAv 1 AQ 1 ) 



AQr, 




Avi-:; 



Figure 5: The relation between the velocity element vfAvi AQq and 
the velocity element u 2 dudVL c . 
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By examining the situation shown in Fig. 5, which is drawn in the velocity 
space, the following relation can be found out: 

/ u 2 dudtt c ■ ■ ■ « vfAvi Afi • • • . (15) 

JAviAflo 

Therefore, the distribution function due to collisions is equal to 

fl 1) (t,r,v 1 ,AQ 1 ) = — \- I dv' I dc' I d{l u , 

2\\J\\a c (u',u)fi°\c> -u>)f!?\c' + u>) (16) 
u\r — r'| 2 

in which u is determined by u = c — v x (vi is in principle along the direction 
of r — r') and u' by u and f2 u /- It should be noted that in (16) differs 
from that in (10) in the sense that AS and Av 1 cease to be arguments and 
V\ takes the place of v x . 

If the zeroth-order distribution functions depend on time and space, the 
replacement 

/fV - u')/ 2 ( V + u') - fi°\t>, r', c' - uU (( V , r', c> + u') (17) 

needs to be taken, where t' — t — |r — r'\/vi stands for the time delay. 

Since AQi has been set finite, expression (16) is nothing but the dis- 
tribution function averaged over A£7 1; which is not good enough according 
to the standard theory. At first glance, if we let AQi become smaller and 
smaller, expression (16) will finally represent the true distribution function 
there. However, the discussion after (12) has shown that if Af^ shrinks 
to zero the definition of the effective cone will collapse and thus the entire 
formalism will no longer be valid. 

A careful inspection tells us that, if we assumed AS* to be finite and 
AQ\ to be infinitesimal at the starting point, a different average distribution 
function, averaged over AS, would be obtained. The fact that no distribution 
function can be determined if Afii and AS simultaneously take on their 
infinitesimal values suggests that even in classical statistical mechanics there 
also exists an uncertainty principle. The connection between this uncertainty 
principle and the uncertainty principle in quantum mechanics remains to be 
seen. 
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In view of that the integration in (16) is carried out over an effective region 
defined by free trajectories of particles, this methodology has been named 
as a path-integral approach[3]. As (16) has partially shown, anything taking 
place in the effective region can make direct impact along 'free' paths and 
any macroscopic structures, continuous or not, will help shape microscopic 
structures elsewhere along 'free' paths. Concepts like the aforementioned, 
though appearing to be foreign from the viewpoint of differential approach, 
are in harmony with what takes place in realistic gases. 

6 Summary 

It has been shown that, when we concern ourselves with beam-to-beam col- 
lisions, part of the system's information has, knowingly or not, been disre- 
garded, and the information loss is characterized by the fact that we have 
four and only four constraint equations (rather than six). Since statistical 
mechanics deals with beam-to-beam collisions, the time-reversibility related 
to each individual collision becomes irrelevant from the very beginning and 
will not reemerge at any later stage. 

By proposing a new context in which what can be really measured in an 
experiment is of central interest, beam-to-beam collisions have been reformu- 
lated. The new formulation reveals that only distribution functions averaged 
over certain finite ranges make good physical sense. 

Apparently, this paper raises many difficult questions related to the very 
foundation of statistical physics. Reference papers can be found in the regular 
and e-print literature [2, 3, 4]. 

Acknowledgments: The author is very grateful to professor Oliver Penrose, 
who reminded me of possibility of defining the cross section of delta-function 
type. Such definition and its possible consequences are discussed in Appendix 
A. This paper is supported by School of Science, BUAA, PRC. 

Appendix A. An alternative cross section 

As revealed in Sect. 3, the original cross section <t(v 1? v 2 — > v[, v' 2 ) defined by 
(2) cannot be interpreted as a function of the usual kind due to the existence 
of the energy-momentum conservation laws. However, one may still wish to 
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define cross sections in which the time reversibility related to a collision of 
two particles plays a certain role. Interestingly, this can be done and doing 
so will help us to find out what kind of problems Boltzmann's equation really 
has. 

If we adopt the definition of (2) and, at the same time, take the energy- 
momentum conservation laws into account, we are led to a cross section of 
delta-function type 

(t(vi, v 2 -> vi, v 2 ) = 77(vi, v 2 -> vi, v 2 )- 

* (VK) 2 + W2) 2 - yfW + M 2 ) ■ 6* [(vi + v 2 ) - ( Vl + v 2 )] . (18) 

To work out the physical meaning of 77, we integrate the above expression 
over a finite but small volume element AviAv 2 



77(vi,v 2 -> vi,Va) 

Av Av' 



* (/WW - v^i) 2 + {v 2 y) ■ s 3 [K + v 2 ) - ( Vl + v 2 )] rfv;rfv 2 . 

(19) 

As shown in Sect. 3, these delta-functions define an accessible shell, denoted 
by S there, and the integration becomes 

77 ( Vl ,v 2 ->vi,v^)AS', (20) 

where AS is enclosed by Av^Av^ (actually by one of Av^ and Av 2 since the 
two are not independent of each other). By comparing this with the cross 
section defined in the center-of-mass system cr c (u, u')d£l, it is clear that 

?7(vi, v 2 -> vi, v 2 ) = <7 c (u, u')/u 2 . (21) 

The above formalism has illustrated that the introduced delta-functions make 
good sense as long as they are integrated over an adequate volume (finite or 
infinitely large). Similarly, 

°"( v i. v 2 -»• v i> v 2) = v(v'i, v 2 -> vi, v 2 )- 

5 ( VK) 2 + W) 2 - VM) 2 + W2) 2 ) ■ 5 3 [(V! + v 2 ) - (vi + v 2 )] . (22) 

With help of (21), we see that 

77(vi, v 2 -> vi, v 2 ) = ?7(vi, v 2 -> vi, v 2 ). (23) 
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This reflects the fact that all collisions, including the original collision and 
the inverse collision, are of head-on type in the center-of-mass system. 

It is now in order to find out whether or not the newly defined cross section 
is relevant to Boltzmann's equation. We first examine the particles leaving 
drdvi during dt because of collisions. Following the textbook treatment, the 
number of collisions is represented by 

[2uf(v 1 )dv 1 ][f(v2)drdv 2 ]<T(y 1 , v 2 -> v[y 2 )dt. (24) 

Integrating it over v' l5 v 2 and v 2 yields, by virtue of (18), (20) and (21), 

dtdvdvxj 2M/(v 1 )/(v 2 ) ( j c (u,u / )rfv 2 rffi. (25) 

Dividing it by dtdrdvi, we obtain the collision number per unit time and 
unit phase volume 

J 2«/(v 1 )/(v 2 )<7 c (u, u')dv 2 da (26) 

If we adopt the assumption that the above number is identical to the number 
of particles leaving the unit phase volume per unit time because of collisions 
(though a different conclusion is offered in Ref. 2), we find that the above 
derivation is entirely consistent with that in the standard approach. 

Then, we examine the particles entering drAvi during dt because of col- 
lisions (Avi has been set finite for a reason that will be clarified). To make 
our discussion a bit simpler, it is assumed that there are only two incident 
beams 

/(vl)Avi and /(v 2 )Av 2 . (27) 

Again, following the standard approach, we know that the collision number 
caused by the two beams within dr during dt is 

[2«/(vi) Avi] • [/K)Mvy • <j(vi, v 2 - Vl , v 2 )dt. (28) 

At this point, a sharp question arises. Can we identify these particles as 
those entering rfrAvi during dt? The answer is apparently a negative one. 
As one thing, only a small fraction of the scattering particles will enter Av 1; 
and the following integration needs to be done: 

/ 2w/(vi) Avi • /(v 2 )drAv 2 • cr(vi, v 2 -> Vi, v 2 )rftrf Vl rfv 2 , (29) 

J Avi Av2 
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in which Avi has to be finite, since the integrand contains delta-functions, 
while Av2 can be infinitely large. If we are interested in knowing the number 
of the entering particles per unit phase volume and unit time, we are supposed 
to evaluate 

AN 

lim - , A , (30) 

where AN is nothing but expression (29). Unfortunately, there are problems. 
As one thing, it has just been pointed out that Avi has to be finite. As 
another, expression (4) and the accompanying discussion have shown that 
the limit expressed by (30) does not exist. 

The above discussion, verifiable with numerical computation, suggests 
that (i) The methodology of determining particles entering a phase volume 
element must be very different from that of determining particles leaving 
a phase volume element, namely there is no symmetry, (ii) Taking limits 
dt — > 0, dr — > and dv\ — > 0, simultaneously or not, can create unexpected 
problems, (iii) In order to formulate the collisional dynamics, new concepts 
and new methodologies are in need. 
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